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ABSTRACT 

We present deep Hubble Space Telescope (HST) NICMOS 2 F160W band observations of the central 
56"x57" (14pcxl4.25pc) region around R136 in the starburst cluster 30 Dor (NGC 2070) located in 
the Large Magellanic Cloud. Our aim is to derive the stellar Initial Mass Function (IMF) down to 
~1 Mq in order to test whether the IMF in a massive metal-poor cluster is similar to that observed 
in nearby young clusters and the field in our Galaxy. We estimate the mean age of the cluster to be 
3 Myr by combining our F160W photometry with previously obtained HST WFPC2 optical F555W 
and F8I4 W band photometry and comparing the stellar locus in the color-magnitude diagram with 
main sequence and pre-main sequence isochrones. The color-magnitude diagrams show the presence 
of differential extinction and possibly an age spread of a few Myr. We convert the magnitudes into 
masses adopting both a single mean age of 3 Myr isochrone and a constant star formation history 
from 2 to 4 Myr. We derive the IMF after correcting for incompleteness due to crowding. The faintest 
stars detected have a mass of 0.5 Mq and the data are more than 50% complete outside a radius of 
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5 pc down to a mass limit of 1.1 M0 for 3 Myr old objects. We find an IMF of 

over the mass range 1.1-20 M© only slightly shallower than a Salpeter IMF. In particular, we find 
no strong evidence for a flattening of the IMF down to 1.1 Mq at a distance of 5 pc from the center, 
in contrast to a flattening at 2 Mq at a radius of 2 pc, reported in a previous optical HST study, 
flattening at 2 Mq at a radius of 2 pc previously found. We examine several possible reasons for 
the different results including the possible presence of mass segregation and the effects of differential 
extinction, particularly for the pre-main sequence sources. If the IMF determined here applies to 
the whole cluster, the cluster would be massive enough to remain bound and evolve into a relatively 
low-mass globular cluster. 

Subject headings: stars: mass function — stars: pre-main sequence — stars: formation — globular 
clusters and associations: individual(30 Dor) 



1. INTRODUCTION 
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The shape of the stellar Initial Mass Function (IMF) 
and whether it is universal or not are key issues in as- 
trophysics. For clusters within 2 kpc, there is no com- 
pelling evidence for variations in the stellar IMF (e.g. 
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IMever et alJ[2000l : IKroupal 120021: IChabried 120051) or the 
brown dwarf IMF (e.g. Andersen et al.ll2008D . However, 
these clusters only span a limited range in total cluster 
mass (10^ — 10^ ^q) and all have a metallicity simi- 
lar to the solar value. Thus, we are forced to observe 
more extreme regions of star formation in search of vari- 
ations in the IMF as a function of environment. It has 
been suggested that the shape of the IMF and in par- 
ticular the characteristic mass where the IMF flattens 
from a Salpeter power law could depend on the metal- 
licity in the molecular cloud out^of w hich the sta rs are 
formed. ILow fc Lvnden-Belll (197^ . rLairemll (|1998f ). and 
lOmukail ()2000[ ) suggest that a lower metallicity results in 
higher temperatures in the molecular cloud which would 
increase the Jeans mass. This would in turn result in a 
top heavy IMF relative to the solar metallicity IMF. 

The closest place with massive metal-poor young star 
clusters is the Large Magellanic Cloud (LMC). The 
metallicity is only | — ^ the solar value (jSmith I [19991 ) 
and star clusters can be studied in some detail despite 
a distance of ~50 kpc (Wcstcrlund 1997). Of particu- 
lar interest is the 30 Dor cluster which is powering the 
most luminous HII region in the Local Group (Kennicutt 
I1984D . The cluster has a i nass of at least 2.2 x10"^ Mq 
within a radius of 4.7 pc ([Hunter et al.|[T995f ) and is a 
relatively low-mass analog to the more distant starburst 
clusters. R136 lies at the center of the 30 Dor clus- 
ter and has long commanded significant attention: Once 
thoug ht to be a single ~1000 M0 star (jCassinelli et al.l 
119811) the region is n ow kn own to host numerous O stars 
(iMelnicklfrgsa IWg elt & Baia1[T985t iPehlemann eTall 
119921 : ICamobell et al.\\im2) . 

The whole 30 Dor region, with a s ize of 200 pc, appears 
to ha v e an age sprea, d of ^20 Myr ([McGregor fc Hvlandl 
119811: iSelman et all [l999i) with stars still forming 
(|Rubio et al.l 119921 : iMaercker fc BurtonI [2005h . R136 
appears to have a much smalle r age spread of at 
most a few Myr (iMelnickl 119851 : iBrandl et aP 119961 : 
iMassev fc Hunted 11998^ An age of 2 Myr or less is 
inferred from spectroscopy of the O stars in the very 
cluster center ([Massev fc Hunteilll998D . whereas the in- 
teri nediate mass popula tion is thought to be ~3-4 Myr 
ol d (iHunter et al.lll995D. 

IMassev fc Hunted (|1998f ) obtained HST spectroscopy 
of the 65 bluest and most luminous sources within 17" 
of the cluster center. They derived the IMF over the 
mass range 15-120 Mq and found it to be well approx- 
imated by a power-law ^fj^^ oc AI^ with a slope of 
F = —1.3 lb 0. 1, c onsistent with a Salpeter slope IMF 
(jSalpeted 119550 . IHunter et all (|1995l Il996f ) obtained 
F555W {V) and F8I4W {i) band optical photometry 
utilizing HST/WFPC2 in order to resolve the cluster's 
intermediate mass stellar population. The IMF derived 
for different annuli out to a radius of 4.7 pc was found to 
be in the range —1.46 < F < —1.17 for the mass range 
2.8-15 MfT), again consistent with a Salpeter slope IMF. 
IMassev fc Hunted (|1998f ) combined their results for the 
high- mass IMF with the results from lHunter et aP ()1995L 
Il996f ) in order to constrain the IMF from 2.8 Mq up to 
120 Mq. Comparing the number of high- mass stars pre- 
dicted by the intermediate-mass IMF from lHunter et ahl 
(|1996[ ). they found the number of massive stars was con- 
sistent with a single power-law IMF with a Salpeter 



slope, i.e. F = —1.35. 

Combinin g the two data sets used in IHunter et al.l 
(fT99l fT99l) . ISirianni et all ((2000h derived the IMF be- 
tween 1.35 Mq and 6.5 Mq, extending the IMF deter- 
mination into the mass range where the stars are still in 
their pre-main sequence phase. The IMF was derived in 
a box with the dimensions ~30'.'4x26'.'8" (7.6pcx6.7pc), 
but excluding the inner most 13'.'6x8.6" (3.5pcx2.2pc). 
Again, a Salpeter slope was found down to 2 Mq, but 
the IMF was found to be flatter than Salpeter, F = 
-0.27 ± 0.08, between 1.35 Mq and 2 Mq, suggesting 
the characteristic mass is higher in this massive, metal- 
poor cluster th an ~ 0.5 Mq as found in the Galactic field 
(|Kroupall200l . 

The foreground (Ay = 0.7 mag) and differential 

extinction (A y 2: — 2 mag) within the cluster 

(I Brandl et al.l [1996) makes it desirable to observe the 
cluster in the infrared, for example the H band where the 
extinction is less than 20% that of the V band. In addi- 
tion, pre-main sequence stars are often associated with 
circumstellar disks and outflows which will introduce ad- 
ditional extinction for the clusters low-mass content. 

We have observed R136 with HST/NICMOS Camera 2 
through the F160W band, which is similar to a ground- 
based H filter. The observations were aimed at being 
sensitive to objects below 1 Mq for a stellar population 
with an age of 3 Myr. Preliminary resu lts have pre- 
vio usly been presen t ed in IZinnecker et al.l ([1999. 2002), 
and lAndersen et all ([20051 ). 

The paper is structured as follows. The data and their 
reduction is described in Section 2. Section 3 shows the 
results for the F160W band imaging. The IMF is de- 
riv ed in Section 4 and compared with the IMF derived 
by ISirianni et al.l ([2000f ) . We point out several plausi- 
ble reasons for the different results in the optical and 
near-infrared, including mass segregation, and differen- 
tial extinction. Finally, our conclusions are presented in 
Section 5. 

2. DATA REDUCTION AND PHOTOMETRY 
2.1. Observations 

We have obtained HST/NICMOS Camera 2 images 
through the F160W band of the central 56"x57" re- 
gion around R136 in the 30 Dor cluster (HST program 
ID 7370). The observations were centered on the clus- 
ter (RA,DEC)=(05:38:43.3, -69:06:08) and on two ad- 
jacent control fields centered on (05:38:42.4,-68:52:00), 
and (05:38:56.9,-68:52:00). The observing dates were 
Oct 14 and 16, 1997. The field-of-view of the 256x256 
pixel NICMOS Camera 2 is 19"xl9" with a pixel scale 
of 0'.'075, resulting in Nyquist sampling of diffraction- 
limited F160W band data. Each position in a 3x3 mo- 
saic centered on R136 was observed four times with small 
dithers of '^16 pixels. The data were obtained in non- 
destructive MULTIACCUM mode such that the photom- 
etry of the bright stars can be retrieved due to the first 
short integration in each exposure. The integration time 
for each dither position was 896 seconds, resulting in a 
total integration time of 3584 seconds for each position 
in the mosaic. The two control fields were observed in a 
similar manner. 

The location of the mosaic is shown in Fig. [1] and the 
NICMOS mosaic is shown in Fig. [2| The faintest stars 
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visible with the stretch used here have an F160W magni- 
tude of ^21.5 mag, corresponding to a mass of 0.8 Mfr^ 
based on the pre-main sequence models of iSiess et alj 
(|2000f ). adopting an age of 3 Myr (jHunter et al.lll995n . 
half solar metallicity, and an extinction of Ay = 1-85 
mag (see Section 3.2). For comparison, the similar de- 
tection limit in an uncrowded environment without neb- 
ulosity would be ~23.5 mag according to the NICMOS 
exposure time calculator. 

2.2. Data reduction 

Each individual image was processed through the 
calnica and calnicb pipelines as well as the biaseq 
and pedsky procedures within the IRAF environment. 
The tasks are described in detail in the NICMOS Data 
Handbook. We used synthetic dark frames and flat fields 
created for the appropriate instrument temperature at 
each exposure. The biaseq task corrects differences in 
bias levels for each chip between different sub-exposures. 
The pedsky task corrects differences in the bias level for 
each quadrant of the chip when the array is reset before 
the exposure. The data for each position in the mosaic 
were combined using the drizzle task. The reduced 
pixel size (0'.'G375) was chosen as half the detector pixel 
size. Bad pixels, bad columns, and the coronagraphic 
hole were flagged as bad pixels before the images were 
combined. 

2.3. Source detection and photometry 

Source detection was done using daofind and pho- 
tometry was performed via point spread function (PSF) 
photometry utilizing allstar within the IRAF environ- 
ment. It was difficult to obtain a good PSF model from 
the data due to the high degree of crowding and the 
spatial v ariability of the PSF . Instead, the TINYTIM 
software ()Hook fc: Kristl[l997l ) was used to create a syn- 
thetic PSF. TINYTIM allows to create a PSF that varies 
as a function of the location on the array. The source 
detection and photometry was performed on each indi- 
vidual position in the mosaic due to the linearly varying 
Point Spread Function (PSF). A hot template star (05V) 
was used for the spectral energy distribution in order to 
achieve the best fit for the brightest stars to limit their 
residuals. The TINYTIM PSF was created for five dif- 
ferent positions on the NICMOS Camera 2 array and 
the PSFs were placed in an empty frame with the same 
number of pixels as the NICMOS Camera 2 array. Four 
frames were created with offsets between each PSF iden- 
tical to the offsets used for the science data in order to 
replicate the data as closely as possible. The four PSF 
frames were then combined using drizzle together in the 
same manner as the science data and a linearly-varying 
PSF was created from the drizzled frame. 

Source detection is complicated due to the diffraction 
features present in NICMOS data. Adoption of a low 
threshold for source detection led to numerous diffrac- 
tion spots from bright stars being erroneously identified 
as fainter stars. Instead the source detection was done in 
the following way in order to limit false detections. We 
first detected the brightest stars (brighter than lOOOcr) 
in each frame and used allstar to remove these with 
the synthetic PSF. A search for fainter stars (brighter 
than 500(t) was then performed in the frame with the 



bright stars removed. Since the removal of the bright- 
est stars also removed the diffraction pattern associated 
with them, we did not detect the diffraction spots as 
stars. The two star lists (the brightest stars and the 
fainter stars) were joined into one and these stars were 
removed from the original frame, again using allstar. 
Fainter stars are then found from the frame with the al- 
ready detected stars removed. This process was iterated 
until stars at 10 a peak pixel intensity over the back- 
ground were detected and removed. The frame with the 
stars removed was then ring-median filtered to remove 
stellar residuals but to retain the large-scale nebulosity 
in each frame. The median-filtered image was then re- 
moved from the original frame and the star detection 
process was repeated in this frame but now continued 
to a detection threshold of 5cr. A ha instead of e.g. a 
3(7 threshold was selected to limit the risks of false de- 
tections due to noise spikes. We finally made sure by 
visual inspection that every detection indeed was a point 
source and that it was not a spurious detection due to 
the diffraction spikes and spots from bright stars. 

The main interest here is in the low-mass (faint) stellar 
content in R136 and one concern is the detection of resid- 
uals from the bright stars as false stellar objects. Some 
false sources were detected by daofind but are rejected 
during the PSF fitting routine. A few remained from the 
brightest stars. They typically produced at most a few 
false detections in the diffraction spikes that were ^--^6-7 
mag fainter than the bright source. We removed these 
detections together with other false detection through 
the visual inspection of all sources. We have further uti- 
lized the artificial star experiments described below to 
examine how many detections are false due to the resid- 
uals from bright stars. We had only false positives as- 
sociated with the brightest artificial stars {F160W < 12 
mag). For artificial stars fainter than F160W^ 14 mag, 
no false detections were present. The false detections for 
the bright stars were located at the diffraction spikes and 
would have been identified in the manual inspection of 
the source list. 

We found a total of 10108 uniquely detected sources 
with a formal error smaller than 0.1 mag and brighter 
than F160W=22.5 mag in the 9 frames. Below this 
magnitude limit the incompleteness is substantial, as dis- 
cussed below. Table [2] presents the list of detected stars. 

2.4. Completeness corrections 

The effects of crowding were examined by placing ar- 
tificial stars in the individual frames using the PSF cre- 
ated from the synthetic TINYTIM PSF. The artificial 
stars followed a luminosity function with a similar slope 
to that of the observed stars (see Section 3) but with a 
surface density 10% that of the detected number of stars 
to avoid affecting the crowding characteristics of the real 
stars. We performed 100 artificial star experiments for 
each frame, for a total of 10 times more artificial stars 
than real stars. Fig. [3] shows the resulting recovery frac- 
tions as a function of the input magnitudes for several 
annuli around the cluster center. The difference of the 
size of the error bars as a function of distance from the 
cluster center is due to a lower number of artificial stars 
placed in the central parts of the cluster. This is a conse- 
quence of adding 10% artificial stars relative to observed 
stars in each artificial star experiment and the relative 
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number of stars in each annulus. The IMF is not de- 
termined in regions with this low completeness. We are 
mainly interested in the low-mass stellar content of the 
cluster, which is below the 50% completeness in the cen- 
tral parts. The uncertainty in the completeness correc- 
tions for the inner parts of the cluster will therefore not 
affect the conclusions drawn for the stellar populations 
further out. 

The completeness is a strong function of the radial dis- 
tance from the center. For the outer regions of the clus- 
ter, 50% or more of the stars brighter than F160W— 21.5 
mag are detected, whereas only the very brightest stars 
are detected in the innermost region. In an annulus at 
0.6-1 pc radius from the center, we detect 50% or more 
of the stars brighter than F160W= 18.0 mag. Adopt- 
ing the PMS models_of Sicss ct al. (2000) and the main 
sequence models of iMarigo et al.l (12008), F 1 60 W= 21.5 
mag corresponds to a 0.8 Mq, half solar metallicity, 3 
Myr old object, whereas F160W— 18 mag corresponds 
to a 7.5 Mq star, assuming an extinction of Ay — 1-85 
mag in both cases. 

2.5. Photometric accuracy 

We have investigated the accuracy of the derived pho- 
tometry by using the stars detected in the overlap regions 
of several fields. Fig. |4] shows the difference in derived 
magnitude for stars detected in the overlap regions of the 
mosaic. Dots denote stars outside a 2 pc radius and plus 
signs denote stars between 1.25-2 pc radius, respectively. 

The F160W band photometry has been compared 
with the ground-based H band pho tometry obtained us - 
ing adaptive optics observations bv lBrandl et al.l (|1996f ). 
Fig. m shows the magnitude difference between the adap- 
tive optics photometry and this study based on 829 stars 
common to both datasets. Stars were considered de- 
tected in both datasets if the spatial position coincided 
within 2.5 drizzled pixels, corresponding to 0'.'094. Some 
scatter is present between the two datasets, especially 
for the fainter stars. However, the median difference be- 
tween the magnitudes derived for the two datasets is less 
than 6% for objects F160W< 18 mag. We have in the 
following treated the F160W observations as a standard 
Cousins H band. 

Conversely, there appears to be a tendency for the 
fainter stars to be b righter in the F160 W data than in 
the H band data of iBrandl erall (11996'). The tendency 
for the fainter stars to be skewed towards fainter H band 
magnitudes is an effect also seen in other comparisons 
betw een HST/NICMOS and AO data (e.g. IStolte eFaH 
120021 ) who suggest it is due to the extended halos present 
in AO observations around bright stars. 

3. RESULTS 

The immediate results from the F160W band HST 
photometry are presented. After discussing the luminos- 
ity function for different annuli, the luminosity profile for 
the cluster is derived. The F160W band data are coni - 
bined with the optical HST data bv lHunter et al.l (|1995l ) 
and the color-magnitude diagrams are presented. Uti- 
lizing the two color-magnitude diagrams we show that 
the spread observed for the higher mass stars is consis- 
tent with that expected due to reddening. We estimate 
the average age for the stellar population and discuss the 
possible presence of an age spread. 



3.1. Luminosity functions 

The star counts in the central 0.6 pc radius region is 
heavily affected by low number statistics, crowding even 
for the brightest stars, and relatively uncertain incom- 
pleteness corrections. We therefore focus on the sources 
outside 0.6 pc in this paper. Fig. [5] shows the F160W 
band luminosity functions for the 0.6-7 pc radius region 
of the 30 Dor cluster divided into several radial bins to 
show the difference in photometric depth due to crowd- 
ing. Overplotted are the completeness-corrected LFs, 
where each bin has been divided by the corresponding 
recovery fraction from the artificial star experiments. 

The completeness-corrected luminosity functions are 
relatively smooth and have been fitted with power-laws 
down to the 50% completeness limit. The derived slopes 
with their la uncertainties and the 50% completeness 
limits are presented in Table [1] Although the slope in 
the inner annulus is found to be more shallow, the de- 
rived slopes are consistent with each other within 2cr with 
an average slope of 0.31 and the shallow slope is not sig- 
nificant. 

The completeness-corrected combined histogram for 
the stars detected in the two off-cluster control fields 
is shown in the lower right panel in Fig. [5l From the his- 
togram, it can be seen that the field star contamination 
found from the star counts is RilO-15% for the faintest 
stars in the 5-7 pc annulus and less closer to the center 
as well as for brighter stars. Although the contamination 
of field stars is found to be relatively small it is not negli- 
gible and they are therefore statistically subtracted from 
the cluster population in the following analysis (Section 
4). 

3.2. The optical-near infrared color-magnitude 
diagrams 

Next, the F160W band phot ometry is combined with 
the optical data presented by iHunter et al.l (|1995D . A 
star was considered detected in both surveys if the spatial 
position agreed within 2.5 drizzled NICMOS Camera 2 
pixels (0'.'094). In the cases where two optical stars were 
located within the search radius of the star detected in 
the NICMOS Camera 2 observations, the brightest star 
was chosen as the match. 

We find in total 2680 in common with the lHunter et al.l 
survey that detected 3623 stars the inner 35" of 
the cluster. 1848 of those sources have a combined formal 
photometric error in the F555W-F160W color of less 
than 0.1 mag. Within the area covered bv lHunter et al.l 
(fT995h we detect a total of 5095 sources. Most of the stars 
detected by the NICMOS survey but not the WFPC2 
observations are fainter than F160W= 20 mag. Assum- 
ing an object age of 3 Myr and an average extinction 
of Ay = 1.85 mag (see below), the similar object would 
have a magnitude in the F8I4W band of ~22 mag. For 
objects with more extinction, they will be even ha rder to 
detect in the F8I4 W band. IHunter et all ()1995D essen- 
tially don't detect any stars within a 1 pc radius at this 
magnitude or fainter. Only 1 in 4 stars in the magnitude 
interval F8I4W =21-22 mag was detected outside 1 pc. 
It is thus not surprising that a significant population of 
faint stars are detected in the NICMOS survey relative 
to the WFPC2 survey. Nevertheless, the lower spatial 
resolution of the NICMOS observations results in a low 
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recovery fraction at these magnitudes in the central few 
pc of the cluster. 

The majority of the sources not detected in the NIC- 
MOS survey but at optical wavelengths are located 
within a radius of 1 pc. The lack of detection is due 
to the lower spatial resolution in this study relative to 
the optical HST data. The resolution is almost a factor 
of two better in the F8I4W band than in the F160W 
band. Sources not detected in the NICMOS data out- 
side 1 pc are mainly due to crowding as well. Indeed, 
visual inspection of the location of the stars detected in 
the optical but not near-infrared shows they are often 
located either very close to the core or on the first Airy 
ring of a bright source. 

The F555W-F160W versus F160W color-magnitude 
diagram is shown in Fig. [6l Overplotted are a 3 
Myr isochrone for the high-mass stars adopted from 
the Marigo et al. [2008.) models and 2, 3 , and 4 Myr 
isochrones adopted from lSiess et al.l (|2000[ ) for stars be- 
low 7 Mq . The stars above 7 Mq and up to the maximum 
mass we fit the IMF in Section 4.1 (20 Mq) are all ex- 
pected to be on the main sequence. Both isochrones were 
calculated adopting a metalli city of half the solar value, 
typical for the LMC (jSmith lil999f ). The two isochrones 
have a small offset in both the V (0. 06 mag) and H (0.07 
mag) band. We have forced the [Marigo et al.l ()2008D 
isochrone to match the iSiess et al.l (|2000f) isochrone at 
7M0. ^ 

It is evident there is a significant scatter in the color- 
magnitude diagram. The scatter is likely due to a combi- 
nation of binary systems (both physical and chance align- 
ments), differential extinction, photometric errors and a 
possible age spread. The median extinction is found for 
the main sequence part of the isochrone. For objects 
in the range 7-20 Mq, we find a median extinction of 
Ay — 1.85 n rag which is sl i ghtly higher than the redden- 
ing found bv lSelman et al.l (|1999f ) in the inner part of the 
30 Dor region. 

At masses below 7 Mq, the spread in the color- 
magnitude diagram is larger but almost exclusively ex- 
tends to the red part of the diagram. This indi- 
cates the lower-mass objects on average have an excess 
amount of ext i nction relative to the higher mass objects. 
iSelman et al.l (|1999( ) observed stars more massive than 
10 Mq and would not detect the additional reddening 
for the lower-mass sources. The possible sources for the 
additional reddening is described in Sec. 3.3. 

We have estimated an average age for the cluster by 
utilising the fact the isochrone is almost horizontal in the 
color range F555W-F160W=1.5-2.5 mag and around 
F160W'^19 mag. The median F160W magnitude is 
19.0 mag in this region of the color-magnitude diagram. 
Adopting an average extinction of Ay — 1-85 mag, this 
corresponds to the F160W magnitudes of the 3 Myr 
isochrone in the same color range. We have thus adopted 
3 Myr as the mean age of the low mass cluster pop- 
ulation and a 3 Myr isochrone is adopted to create a 
mass-luminosity relation in order to turn the luminosi- 
ties into mas s es for objects below 7 M© and the 3 Myr 
iMarigo et al.l (|2008f ) isochrone above. We will in Sect. 4 
the effects on the derived IMF adopting an age spread of 
2 Myr. 

The right hand panel in Fig. [6] shows the I-F160W 
versus F160W color-magnitude diagram. It is evident 



that the clustering around the isochrone is tighter than 
for the V-F160W versus F160W color magnitude dia- 
gram. This is expected if a large part of the scatter is 
due to differential extinction. We can calculate the scat- 
ter around the main sequence in both color-magnitude 
diagrams and compare with the difference predicted from 
extinction. 

If the spread in the color-magnitude diagrams is due 
to extinction we expect the ratio of spread in the V- 
F160W versus F160W diagram to be ratio of the extinc- 
tion in each color, i.e. (1 - 0.192)/(0.62 - 0.192) 1.88 
times larger than in the I-F160W versus F160W color- 
magnitude diagram. Since the isochrone is almost ver- 
tical in both diagrams, we have calculated the standard 
deviation around the reddened isochrone in both color- 
magnitude diagrams. We have used the stars with good 
photometry, better than 5% in each filter, and in the 
magnitude range 13 < F160W < 17mag. The standard 
deviation found for the V-F160W and I-F160W color- 
magnitude diagrams are 0.60 mag, and 0.36 mag and the 
ratio is 1.7. If the measurement errors are taken into 
account this ratio increases. The typical errors for the 
culled sample are 0.04, 0.02, and 0.03 mag for the V, I, 
and F160W bands, respectively. After taking the mea- 
surement errors into account, the ratio is found to be 
1.9, assuming the measurement errors in two filters are 
independent. Due to blending, this is not necessarily the 
case. Thus, 1.9 is an upper limit and we thus find the 
ratio to be between 1.7 and 1.9, in agreement with the 
scatter being due to differential extinction. 

Since the amount of differential extinction does not af- 
fect the F160W band photometry significantly, the sin- 
gle band photometry presented here is competitive with 
the 2-band optical photometry. There is a unique trans- 
lation from the F160W band magnitude to the object 
mass for the majority of the mass range. For the optical 
photometry, the color information is used to determine 
the extinction and the mass function is thus effectively 
determined by the de-reddened V band magnitude. 

3.3. The differences between the optical and 
near-infrared HST observations 

The main advantage of the optical relative to the near- 
infrared HST photometry is the improved resolution due 
to the smaller diffraction limit. The stellar content can 
therefore be resolved to lower masses closer to the clus- 
ter core than is possible with the near-infrared observa- 
tions. However, phenomena associated with the star for- 
mation process can introduce additional reddening that 
can complicate the derivation of the low-mass IMF from 
optical data. The low-mass objects may still be as- 
sociated with a circumstellar disk. There is evidence 
from e.g. the Orion Nebula Cluster that circumstellar 
disks can survive the UV radiation from massive stars 
(jRobberto et al.ll2004[ ) . Even if the disks are being evap- 
orated by the ration field from the early type stars, the 
evaporated material will be a further source of reddening. 
Patchy extinction associated with the 30 Dor complex 
and located in the foreground of R 136 will be an addi- 
tional source of differential reddening. Th ere are signs in 
the op tical images presented in Fig. 1 of iSirianni et al.l 
(|2000r ) of patches of extinction, e.g. to the east-north- 
east of the cluster center. If variable extinction is present 
or if a significant fraction of the stars are associated with 
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disks or outflows, an extinction limited sample has to be 
created in order to avoid a biases against detection of the 
low-mass stars. 

The near-infrared photometry is effected by differen- 
tial extinction as well but the effect is less than 20% of 
that measured in the V band. Thus, whereas the IMF 
derived from optical observations where an extinction 
limited sample is not defined might be severely affected 
for the low-mass objects, the effect on near-infrared ob- 
servations is modest. Therefore, in the outer parts of 
the cluster where crowding is a smaller issue than closer 
to the center, the near-infrared observations are more 
suitable to detect and characterise the low-mass stellar 
population in the cluster. 

On the other hand, single band photometry has the 
disadvantage that there is no information on the age of 
individual objects. We investigate in the next Section 
how this might affect the derived IMF. We note that if 
differential extinction is present, the situation is no better 
for the optical photometry. Even though the cluster was 
observed through two filters in the optica l, there is still 
a degeneracy between age and extinction. iSirianni et al.l 
(jlOOO) converted the V-i photometry into an effective 
temperature and used that effective temperature to ob- 
tain a bolometric correction. Without de-reddening the 
sources, the age of a cluster member can be in error and 
hence the mass estimates will be uncertain. 

4. ANALYSIS 

We construct a mass-luminosity relation by combin- 
ing the main sequence models by Marigo ct al. (2008), 
and the pre-main sequence models of Siess ct al.. (2000,) 
in order to infer the stellar mass from the F160W band 
magnitude. We then derive the mass functions for R136 
outside 0.6 pc where the 50% completeness limit corre- 
sponds to a stellar mass below 10 M©. Der iving the IMF 
this way is a well establi shed procedure (|Lada fc Ladal 
120031: IMuench et al.|[200l ). We further discuss the po- 
tential effect of extinction on the derived IMF. Finally, 
we search for evidence for mass segregation in the outer 
parts of the cluster using the cumulative luminosity func- 
tions. 

4.1. Deriving the mass function 

A mass-luminosity relation is needed to convert the 
derived F160 W band magnit ude for each star to a mass. 
We use the Si ess et al.l (120001 ) isoch rones for stars below 
7 M© and the IMarigo et al.l (|2008l ) 3 Myr isochrone for 
the more massive stars as discussed in Section 3.2. The 
age of the cluster is first assumed to be 3 Myr and is 
later varied to examine the effects on the derived IMF 
for different cluster ages. Stars below '^3 Mq are on the 
pre-main sequence isochrone whereas the more massive 
stars up to our upper mass limit of 20 Mq (see below) 
are on the main sequence. The adopted mass-luminosity 
relation is shown in Fig. [T] 

We have limited knowledge of the extinction for the 
majority of our objects. Instead, we have adopted an 
average extinction of Ay = 1.85 mag, as determined from 
the V-H versus H color-magnitude diagram in Fig. [51 
Since the am ount of extinction ranges between Ay — 
0.7 — 3 mag ([Brandl et al.l [l996l ) the extinction for an 
individual object might be wrong by up to Ay ~ 1 mag, 
a maximum error < 0.2 mag in the F160W band. This 



corresponds to an error of '^10% when the luminosity is 
transformed into a mass. 

Fig. [8] shows the derived mass functions outside 0.6 pc 
for a 3 Myr isochrone after field stars have been sub- 
tracted statistically in each annulus. The mass functions 
are in general smooth and well fit by power-laws. How- 
ever, there appears to be some structure in the derived 
IMFs at intermediate masses, 2-4 M©, which is the re- 
gion where the pre-main sequence track joins the main 
sequence. The mass-luminosity relation is plagued by a 
non-monotonous feature at this mass rang e (sec Fig. [7t, 
which marks the radiative-convective gap (jMavne et al.l 
l2007f ) and the transition r egion from pre-ma in sequence 
to main sequence (see also lStolte et al.l [20041 ) . A similar 
structure in the derived IMF is seen in the results from 
e.g. NGC 3603 but at a slightly higher mass since the 
cluster is younger (e.g. IStolte et al.N2006l ). The turn-on 
mass is higher for a younger cluster. Thus we would ex- 
pect the kink in the mass-luminosity relation to move to 
higher masses for a younger cluster. Since this is what is 
seen comparing NGC 3603 and R 136, it indicates indeed 
a feature of the isochrones and not a feature intrinsic to 
the cluster. The number of stars in each mass bin is 
provided in Table [H . 

Power-laws have been fitted to each of the histograms 
in order to derive the slopes of the mass function in each 
annulus. The fit was done over the mass range from 20 
Mq down to the 50% completeness limit for each annulus. 
The mass for stars above ~20 Mq is very poorly con- 
strained from near-infrared observations due to uncer - 
taintics in the bolometric corrections (e.g. lMassevll2003l ) . 
The derived slopes F, where dN/dlogM oc M^, are in- 
dicated in Fig. [51 and are also presented in Table [31 The 
derived slopes for annuli outside 1 pc are consistent with 
each other within 2a error bars. For the 3-5 pc and 5-7 
pc annuli where the data are complete to below 2 Mq, 
the slopes are found to be —1.2 ± 0.1 and —0.9 ± 0.2, 
respectively, slightly shallower than the s lope of F =- 
1.28 ±0.05 derived by ISirianni et al.l (|2000l ) above 2 Mq, 
except that in our case the IMF continues as a power-law 
down to 0.8 Mq. 

Has the fact that we used the whole mass range for 
our power-law fit washed out a possible flattening at the 
low mass end? To test this possibility, we have addition- 
ally fitted a separate power-law to the low-mass part of 
the IMF. Only the part of the mass function that is not 
influenced by the kink in the mass luminosity relation 
is used. This region is limited to masses below 1.7 Mq 
for the 3 Myr isochrone. It is therefore only for the 5-7 
pc annulus that a reasonable mass range is covered to 
fit the IMF. We find the slope to be F = -0.9 ± 0.2, 
which is more shallow but consistent at the 2a level with 
a Salpeter IMF and is consistent with the slope derived 
for the full mass range. 

We have der i ved th e IMF in the same boxes as done by 
ISirianni et al.l (|2000f) . The completeness correction was 
calculated independently for each box before the IMFs 
were combined to the average I MF for direct c ompar- 
ison with the IMF presented by ISirianni et al.l p0003 . 
The 50% completeness limit for the NICMOS data varies 
from 2.8 to 1.4 Mq for the four boxes. Following Sirianni 
et al., we have derived an average completeness limit for 
the three regions of 2.2 Mq. As evident, the agreement is 
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good for the common mass range. We ap pear to under- 
estima te the stars at ~ 6M0 compared to lSirianni et alj 
(|2000f ). However those appear to be recovered at 8 Mq. 

The color-magnitude diagrams show a large spread in 
the main sequence to pre-main sequence transition at 
F160W=18-19 mag. Although shown in Section 3.2 that 
this scatter can be explained by differential extinction, it 
cannot be ruled out that there is an age s pread present 
as well as suggested in prev ious studies (jHunter et al.l 
ll995HMassev fc Huntei1ll998l ). It is therefore reasonable 
to take a star formation history different than a single 
burst at 3 Myr into account. Wc show in Fig. [9] the 
IMF in the outer two annuli assuming a cluster age of 
2 and 4 Myr, respectively. We also show an 'average' 
IMF found as the average of the IMF's derived for the 
age range 2-4 Myr in 0.5 Myr increments. The lower 
mass limit in the average IMF was determined from the 
4 Myr isochrone which provides the most restrictive mass 
limit. We find that both in the case of a 2 and 4 Myr 
isochrone the IMF is well fit by power-laws. The derived 
slopes are steeper assuming an older isochrone relative to 
the younger ones. There is no indication for a flattening 
below 2 M0 in either case. The average IMF is also found 
to be represented by a power-law with a slope consistent 
with a Salpeter slope. The slopes of the derived power- 
laws are given in Table [31 For the average IMF, the 
number of stars averaged over the different ages in each 
mass bin is derived. Error bars for the average IMF have 
been determined as the standard deviation around the 
mean number of objects in each mass bin. 

As was the case for the 3 Myr isochrone, the slopes 
of the IMF for different assumed ages have also been 
calculated and are provided in Table [S] The slopes are 
found to be shallower than a Salpeter slope, but at the 

2a level consistent with a Salpeter slope. The slopes 
are also consistent with those derived for all masses up 
to 20 M0 , as was the case assuming the 3 Myr isochrone. 

The lack of a flattening in the I MF below 2 Mp is i n 
contrast to the results presented bv lSiriarLifi "eFalllBoOl) . 
who derived the IMF closer to the cluster center. There 
can be several possible reasons for the difference in the 
derived IMF slope in the two surveys. First, due to the 
different spatial resolution in the two studies, the NIC- 
MOS IMF is derived further away from the cente r of th e 
cluster than the WFPC2 IMF bv iSirianni et al.l (j2000D . 
The IMF was derived in the areas shown in Fig. [3] as re- 
gions B,C, and D. Thus, all of their surveyed area is out- 
side a radius of 1 pc and the majority of their surveyed 
area is between 2 and 5 pc where crowding precludes 
NICMOS from detecting stars less massive than 2.2 Mq 
for a 3 Myr isochrone. One possibility for the difference 
in the derived slopes for the NICMOS and WFPC2 data 
can therefore be a variation of the IMF as a function 
of radius. Another possibl e reason can be diff erential 
extinction as suggested by iSelman et al.l (|1999l ). Both 
possibilities are discussed in subsections 4.2 and 4.3 

4.2. The effect of differential extinction 

As was suggested by ISelman et al.l (|1999[ ). the pres- 
ence of differential extinction can potentially alter the 
low-mass end of the IMF if an extinction-limited sam- 
ple is not used. In order to estimate the possible effect 
on the IMF if differential extinction is not taken into 
account, we have constructed a simple model of the clus- 



ter which includes d ifferential extin ction and the depth 
of the dataset from ISirianni et al.l (|2P00). A Salpeter 
slope IMF and a cluster age of 3 Myr were assumed. 
Each object within the artiflcial cluster was assigned 
a V band magnitude based on its mass from the 3 
Myr isochrone co mputed for a half solar metallicity by 
iSiess et al.l (|2000f ). The objects were then reddened by 
a foreground extinction chosen randomly from a normal 
distribution with a standard deviation of 0.7 and shifted 
to peak at Av=1.85 mag. If the extinction was found to 
be less than Ay = 0.7 mag, a new extinction was calcu- 
lated. Stars were then considered detected if their red- 
dened magni tude is within the 50 % completeness limit 
presented bv ISirianni et al.l ()2000D . The model is obvi- 
ously an oversimpliflcation of the real situation. Never- 
theless it is expected to illustrate how the derived IMF 
might differ from the underlying IMF. 

Figure [To] shows the input Salpeter IMF (solid line), 
the derived I MF (dashed li n e) tog ether with the mea- 
surements by ISirianni et al.l ()200Gl ). The model mimics 
a fl attening in the observ ed IMF similar to that deduced 
by ISirianni et al.l ()2000l ) . The ratio of the number of 
stars below and above 2 Mq respectively has been cal- 
culated both f or the model cluster and the data from 
ISirianni et al.l (|2000f) . For the model cluster it is found 
to be 0.87, which is in reasonable agreement with the 
ratio of 0.76 derived from the observations. 

4.3. Cumulative mass functions in the outer parts of 

R136 

Another explanation for the differenc e between the 
result s obtained here and the results by ISirianni et al.l 
(|2000D can be mass segregation. We have searched for 
evidence for mass segregation in the two outer annuli in 
our survey. We used the luminosity functions instead of 
the mass functions to avoid additional uncertainties due 
to the mass-luminosity relation. The results obtained 
for the mass functions are very similar to those from the 
luminosity functions. 

The cumulative luminosity distributions are shown in 
Fig. [m for the outer two radial bins. These are the only 
bins where the 50% completeness limit is below 2 Mq. 
The two cumulative distributions are very similar. We 
have performed a Kolmogorov-Smirnov test to quantify 
the similarity of the cumulative luminosity distributions. 
The maximum difference between the two distributions is 
0.039 and the probability for the two distributions to be 
drawn from the same parent distribution is 10%. Thus, 
there is no strong evidence (less than 2a) for mass seg- 
regation in the outer parts of the cluster. 

The fact that there is little evidence for mass segrega- 
tion outside 3 pc does not exclude the possibility that 
the cluster is mass s egregated out to a rad ius of several 
pc. Both iMalumut h & Heap ( 1994) and iBrandl et all 
(1996) found evidence for mass seg regation of the mas - 
sive stars in the center of the cluster. IBrandl et al.l (|1996l ) 
showed the half-mass relaxation time to be 7.8T0'' yr, 
much longer than the cluster age. They also point out 
that the massive stars will experience mass segregation 
on a much shorter time scale than the lower mass stars; 
the time scale depends inversely on the stellar mass. It 
is thus not surprising, from a dynamical point of view, 
that there is no evidence for mass segregation outside the 
half-mass radius of 1.7 pc (^Hunter et al.„1995l . 
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On the other hand, this does not rule out the 
possibiUty that the cluster might be mass segre- 
gated at birth closer to the cluster center. Evidence 
for mass segregation has b een found in e.g. the 
Orion Nebula Cluster (ONC) (Hillenbrand & Hartmann 
[r998: Bo nnell fc Dav ics 1998). Hillenbrand & Hartmann 
(jl998l ) showed evidence for mass segregation down to 
stellar masses of 1-2 Mq. Due to the youth of the ONC, 
they concluded the mass segregation had to be at least 
partly primordial. It is thus possible that R136 is also 
affected by primordial mass segregation close to the clus- 
ter center and that mass segregation is the reason for the 
difference between the NICMOS and WFPC2 IMFs. 

4.4. Cluster mass 

We can obtain a rough estimate of the cluster mass 
from the near-infrared observations. The main limita- 
tion in our mass estimate is the amount of confusion due 
to crowding in the cluster centre: our data mainly sam- 
ples the IMF down to and below 1.4 M0 outside 3 pc. 
Neverth eless, we can utilise the mass estimates within 2 
pc from lHunter et ahl (|T995) to complement our mass es- 
timate down to 2.1 M0. Inside 2 pc and into 0.15 pc, the 
results are extrapolated from the local completeness limit 
mass down to 2.1 M0 assuming an underlying Salpeter 
IMF. No stars have been detected less massive than 20 
Mq within the central 0.15 pc radius due to crowding. 
The mass in the very center has been estim ated from the 
surface density profile down to 2.8 Mq in iHunter et al.l 
(| 19951 ) to be 4 • lO'* Mgpc'^, resulting in a mass of 3700 
Mq down to a lower mass limit of 2.1 Mq. We find the 
cluster total mass down to 2.1 Mq to be 5 • 10'* Mq. The 
directly determined mass down to 2.8 Mq within 4.7 pc 
is found to be 2 . • 10" ^ Mq, almost the same as found 
by IHunter et all (|1995f) . If the IMF follows a Salpeter 
slope down to 0.5 Mq as obser ved in the Ga lactic field 
and nearby lower-mass clusters (|Kroupall2002l ) , the total 
mass in the central region would be roughly double the 
amount given above, and the total cluster mass would be 
close to ~ IO^Mq. 

The velocity dispersion, and hence the dynamical 
mass, of the whole NGC 2070 region , including R 136 
has been determined by iBosch et al.l (|2009f ). The dy- 
namical mass was determined to be 4.5T0^Mq, almost 5 
times higher than expected for R 136 alone, but con- 
sistent with the ph otometric mass for the same area 
(|Selman et al.llT99l . If we take into account that the 
half mass radius of R 136 is 1.7 pc (jHunter et al.lll995l ). 
compared to 14 pc for the whole NGC 2070 region and 
assuming the velocity dispersion is the same in the inner 
parts of the cluster, we would expect a dynamical mass 
of 4.5-10^ • 1.7/14Mq=5.5-10'*Mq which is lower than the 
mass expected if the IMF is consistent with a Galactic 
IMF down to 0.5 Mq. Thus, at face value, the velocity 
dispersion would be low enough that the cluster can stay 
bound. However, a measurement of the velocity disper- 
sion for the inner regions is necessary to directly compare 
the photometric mass with the dynamical mass. 

4.5. The surface brightness profile 

We can directly derive the surface brightness profile 
of the region around R136 in the 30 Dor cluster since 
the data does not suffer from saturated stars. Although 



bright stars will saturate through the one hour expo- 
sure, the non-destructive readout mode ensures that only 
the first reads are used to derive the magnitude of the 
brightest stars. The surface brightness profile is shown 
in Fig. M 

Between ~0.2 and 2 pc, the light profile is well fit by a 
power-law, whereas inside 0.2 pc the light profile appears 
to be flattening. We have therefore fitted the light profile 
with a power-l aw modified by a c ore radius, similar to 
the approach in lElson et al.l (|1987l ). Constraining the fit 
to inside 2 pc, we find a slope of —1.54 ± 0.02, slightly 
more shallow than —1 .72 ± 0.06 derived outside O.lpc by 
iCampbell erahl (|1992f) using F336W Planetary Camera 
onboard HST observations. 

The core radius is found to be 0.025 ± 0.004pc, which 
is less than the resolution of the observations and is thus 
likely and upper limit. Previous HST optic al studies de- 
term ined a small core radius, rc < 0.02 pc (jHunter et al.l 
|1995[) . consistent with our findings here. However, since 
the derived core radius is smaller than the resolution of 
the observations, it's evidence is weak. One or two bright 
stars off center by only a small amount could mimic a 
cluster core. 

4.6. Comparison with other massive clusters and the 
implications of low-mass stars in R136 

How does the low-mass end of the IMF in 30 Dor com- 
pare with that determined for other massive and dense 
stellar clusters? A top-heavy IMF in massive dense clus- 
ters h as been suggested on theoretical grounds (e.g. ISilkl 
|1995[) . The most convincing example of a young cluster 
with a present-day mass function departing significantly 
from a Salpeter IMF above 1 Mfn is the Arches cluster 
(|Stolte et al]l200llFiger eFalllTggl . fStolte et al.l([200l ) 
found an average slope ofF = — 0.9±0.15 for the central 
parsec of the Arches cluster, flatter than a Salpeter slope 
of -1.35. Deeper observations found that the present day 
mass function in Arches to be well approximated by a 
power-la w with a slope o f F = —0.91 ± 0.08 down to 
1.3 Mq (|Kim et al.l 120061 ). However, recent work tak- 
ing differential extinction into account suggest the slope 
of the power-law is only slightly more shallow than a 
Salpeter slope, F = -1.1 ± 0.2 (jEspinoza et al.l 120091) . 
iPortegies Zwart et al.l ()2002D note that even if the ob- 
served IMF is slightly flatter thann a Salpeter IMF, this 
can be explaned by mass segregation. The mass segrega- 
tion would be accelerated in the cluster due to the strong 
gravitational field from the Galactic Center. By adopting 
realistic parameters for a model cluster and an appropri- 
ate distance from the Galactic center, they found that 
an input Salpeter slope IMF would be transformed to 
the observed prese nt day ma ss function via strong dy- 
namical evolution. IStolte e~al.l (2006) showed that the 
IMF of the cluster powering the NGC 3603 HII region 
was well fitted by a power-law but with a slope flatter 
than Salpeter, F = -0.91 ± 0.15. They further showed 
evidence for mass segregation for the more massive stars, 
M> 4 Mq. The data indicated a slight flattening of the 
low-mass content (M< 3 Mq). NGC 3603 is younger 
than the Arches cluster and not affected by a strong 
tidal gravitational field. Thus it is expected to be less 
influenced by dynamical mass segregation. 

The even more massive starburst clusters appear to be 
the primary sites (unit cells) of star formation in star- 
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burst galaxies, including interacting/colliding galaxies 
such as The Antennae or The Cartwheel. If starburst 
clusters are the basic building blocks of certain star- 
forming galaxies, their stellar content (IMF) will affect 
much of the observed chemical and photometric evolu- 
tion of galaxies, both at the present epoc h and perhaps 
even more so in the high-redshift past (jCharlot et al.l 
|1993[) . Several observational claims have been made that 
the IMF in unres olved starburst clusters is top-heavy 
(jRieke et al.lll993l ). although observations of the Anten- 
nae gave a mixed result (Mcngcl ct al. 2002). However, 
it has been suggested that the high mass-to-light ratios 
found in some young starburst clusters are artificially 
high related to their not being in v irial equilibrium due 
to ga s expulsion from the clusters (|Goodwin fc BastianI 
|2006[) . During the first 50 Myr of the cluster, the ve- 
locity dispersion and hence the cluster mass might be 
overestimated if the cluste r is assumed to be virialized. 
iGoodwin k, Bastiaiil (|2006[ ) suggest that the top-heavy 
IMFs inferred in young unresolved extragalactic star 
clusters might be spurious due to their non-virialized dy- 
namical state. 

With the present dataset it is clear that the IMF in 
the outer parts of R136 continues as a power-law down 
to 1 Mq, similar to what is found in other star clusters 
and the slope is similar to what is found in the field. 
Whether this is true for the cluster as a whole depends 
on the cause for the flattening observed closer to the clus- 
ter center. It would be interesting to know the IMF if 
the observations could be extended closer to the charac- 
teristic mass where the Galactic field star IMF flattens 
(0.5 Mq IKroupal[200l) . a mass that can be reached in 
massive young clusters (< 4 Myr) in the LMC with AO 
systems. 

It has long be en suggested R136 inight be a proto- 
globular cluster (|Mevlanill99"5 ILarsonlll993l ). The ques- 
tion has been whether R136 would remain bound over 
a Hubble time. One consequence of a top-heavy IMF 
is that the cluster would dissolve soon after gas expul- 
sion and mass loss due to evolution of the high-mass 
stars. However, the detection of stars in R136 less mas- 
sive than 1 Mq gives the first direct evidence that low- 
mass stars are formed in a starburst cluster. The fact 
that the IMF in the outer parts of R136 appears to be 
a Salpeter IMF down to at least 1 Mq gives support to 
the notion the cluster might be a proto-globular cluster, 
albeit a light one. Early gas expulsion and subsequent 
mass loss through stellar evolution will disrupt star clus- 
ters deficient in lo w-mass stars during the first 5 Gyr 
of the clusters life (jChernoff fc Weinberelll990t iGoodwinI 
Il997f) However, a determination of the velocity disper- 
sion in the inner parts of the cluster is necessary to de- 
termine its final fate. Thus, the presence of low-mass 
stars is a necessary, but not sufficient condition for the 
possibility of the cluster to evolve into a globular cluster. 
The median mass of Galacti c globular clusters is 8.1-10'' 
Mq ()Mandushev et anil991f ). comparable to the mass of 
R136. Even if R 136 will remain bound it will lose some 
mass and might end up as a low-mass globular cluster. 

5. CONCLUSIONS 

We have analyzed HST/NICMOS F160W band data 
covering the central 14pcxl4.25pc around R136 in the 
NGC 2070 cluster in the LMC. We have reached the fol- 



lowing conclusions: 

• From the color-magnitude diagram obtained by 
combining our photometry with previously pub- 
lished HST/WFPC 2 F555W data we constrain 
the age of the lower-mass stellar content in the 
cluster to be 2-4 Myr, consistent with previous es- 
timates. We derive individual masses for the ob- 
jects detected adopting a 3 Myr isochrone. 

• We have detected stars in the cluster down to 0.5 
Mq at r > 5 pc, assuming an age of 3 Myr. 

• The derived IMF is consistent with a Salpeter slope 
IMF with no evidence for a flattening at low masses 
down to the 50% completeness limit corresponding 
to a mass of 1.1 Mq outside a radius of 5 pc for a 3 
Myr population and 1.4 Mq if the oldest stars are 
4 Myr. 

• The result is in disagreement with the flattening 
of the IMF below 2 Mq observed bv lSirianni et al.l 
(|2000h using optical data covering a region closer 
to the cluster center. We suggest two possible rea- 
sons for the discrepancy: differential extinction and 
mass segregation. 

• We find no evidence for mass segregation outside 3 
pc, but with the current data, we cannot rule out 
that closer to the center the low-mass stars are 
segregated. 

• From the radial surface brightness profile we have 
derived a core radius for the cluster of 0.025 
pc (O^^l), consisten t with previous estimates by 
iHunter et all (fl995l) . 

• The mass of the cluster within 7 pc between 25 Mq 
and down to 2.1 Mq is estimated to be 5T0^ Mq. 
If the IMF continues with a Salpeter slope down to 
0.5 Mq the total mass estimate will double. 

• The total mass of the cluster combined with the 
large number of low-mass stars suggests that the 
30 Dor cluster may survive to become a proto- 
globular cluster depending on the cluster velocity 
dispersion. 
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Fig. 1. — Location of the NICMOS surveyed region within 30 D or, outhned by t he black box, overlaid on a, K band image obtained 
with the IRAC-2a camera on the MPI/ESO 2.2 meter telescope (see lBrandner et al.ll200ii . for details). The full field is 200" x225"and the 
NICMOS field is 56"x57". North is up and East is left. 



12 



Andersen et al. 




Fig. 2.— The central region around R136 in the 30 Dor cluster (NGC 2070) as observed through the F160W filter with the NICMOS 
Camera 2 on HST. The field-of-view is 56"x57", corresponding to 14 pc X 14.25 pc. A logarithmic inverted intensity scale has been used. 
North and Ea st are indicate d in the Figure. The faintest sources visible in the image are F160W ~21.5 mag. Stars with spectral types 
determined bv lParlterl II1993I ) are marked by their ID number, along with R134, MK34, and MK42 |Melnick 1985) to the lower right of each 
star. All the identified stars are early O or WR type stars, except P1257 (BOI A), P1253 (BN0.5Ia) , and P987 (B0.5-0.7I). The spectral 
types are adopted from Walborn & Blades (.1997 ). The four regions analysed bv lSirianni et aLl 12000) are shown and labeled. Further, the 
1 pc and 5 pc radii are indicated by circles. 




Fig. 3. — The completeness as a function of magnitude for different annuli extending from the center to the edge of the field. The error 
bars indicated on every second curve are the statistical fluctuations in the completeness tests. The corresponding masses for a 3 Myr 
isochrone and an extinction of Ay = 1.85 mag are indicated at the top of the plot. 
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Fig. 4. — Left: A comparison of the NICMOS 2 photometry derived in this study with the ground-based AO H band observations by 
IBrandl et al.l II1996I ). The mean difference between the two datasets for objects with F160W < 17 mag is 0.03 mag, and the standard 
deviation is 0.2 mag. Right: The magnitude difference for sources detected in two different images in the mosaic. A star was identified in 
both images if its position was within one pixel. Plus signs denote stars within a 2 pc radius of the center, but outside 1.25 pc; the dots 
denote stars outside this annulus. The RMS in the error for objects with F160W< 21.5 mag, which is the 50% completeness limit in the 
outskirts of the cluster, is 0.17 mag. 



The low-mass IMF in the 30 Doradus cluster 15 



0.60 to 1.00 pc 



.00 to 2.00 pc 



o 


1000 





CL 






& 


100 





V) 








10 





Q_ 








1 





O 






-t— ' 
(/) 





1 




10 12 14 16 18 20 22 10 12 14 16 18 20 22 



o 


100.00 r 


CL 




6" 


10.00 r 


W 




I- 
(U 


1.00 r 


Cl 




in 

V- 


0.10 r 


O 
-t— ' 




CO 


0.01 ? 



F160W band magnitudes 
2.00 to 3.00 pc 



F160W band nnagnitudes 
3.00 to 5.00 pc 




10 12 14 16 18 20 22 10 12 14 16 18 20 22 
F160W bond magnitudes F160W band magnitudes 



o 


100.00 r 


Q_ 




& 


10.00 r 


(f) 






1.00 r 


Q_ 






0.10 r 


O 




CO 


0.01 r 



5.00 to 7.00 pc 



LF for the control field 




10 12 14 16 18 20 22 10 12 14 16 18 20 22 
F160W band magnitudes F160W band magnitudes 



Fig. 5. — The F160W—h&i\d luminosity functions for different annuli in R136 outside a radius of 0.6 pc. Solid histograms indicate the 
observed number of stars per magnitude bin while the completeness-corrected data are shown as the dotted lined histograms. The error 
bars include the Poisson errors and the uncertainty from the incompleteness calculations where the two error terms have been added in 
quadrature. In each panel, an arrow indicates the bin where the completeness correction is 50%. The solid straight line is a weighted fit 
to the completeness corrected histograms down to the 50% completeness limit. The lower right panel shows the completeness corrected 
average luminosity function for the two control fields (solid histogram) associated with the 30 Dor cluster. The completeness corrected 
luminosity function for the 5-7 pc annulus in 30 Dor is shown as the dashed histogram for comparison. 
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Fig. 6. — Left: The F555W-F160W versus F160W color-magnitude diagram for R136 for the central 35"x35" observed bv lHunter et al.l 
II1995I) in the F555W and F814 W bands. The arrow illustrates the effect of A y = 1 ma g of extinction (Apmow = 0.175X Ay). Overplotted 
as the solid lines are the 2, 3, and 4 My r pre-main sequence isochrones from lSTess et al. ( 2000) below 7 Mq. The long dashed line is the 
main sequence from Marigo et al.l l)200g ) below 7 Mq and the 3 Myr isochrone from a Marigo et al. (2008) above 7 Mq. All the isochrones 
have been reddened by Ay = 1.85 mag, the median reddening for the high mass stars (sec the text) and shifted by a distance modulus of 
18.5. The F5 55W magnitudes have been transformed into the V magnitude using the transformations in Holtzman ct al. (1995). The 3 
Mvr lSiess et al.. (i200(3i ) isochrone joins the main sequence at F160W ^19.5 mag, V—F160W=0.8 mag. Masses arc indicated on the 2 Myr 
isochrone. Right; Same as to the left but this time the l-F160WveTsus F160W color-magnitude diagram. Notice the tighter clustering 
around the isochrone to the right, indicating part of the scatter in the W-F160W color— magnitude diagram is due to differential reddening. 
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Fig . 7. — The mass-luminosity relation for a 3 Myr isochrone. Marked is the location where the ISiess et al.l II2000I ) and IMarigo et aLl 
1)20081 ) isochrones merge (7 M0). The feature at 2—3 Mq is due to the pre-main sequence to main sequence transition. 
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Fig. 8. — The mass functions for 30 Dor outside 0.6 pc in several annuli after field star subtraction. The arrows indicate the 50% 
completeness limit. Dotted line histograms show the mass functions derived from the uncorrected star counts whereas the solid line 
histograms are the completeness corrected mass function. Expected errors due to Poisson noise are indicated on the solid line histograms. 
A fit is made to the completeness corrected histogram, corrected for field star contamination. The maximum mass used in the fits is 20 
M0. The coefficient shown in each panel is F, dN/dlog M oc . The lower right panel shows the IMF derived bv ISirianni et aLl I I2000I ) 
from the ar eas shown in Fig. ^ Further, the plus symbols show the IMF derived from the NICMOS data from the same regions as was 
used in the Itjirianni et aTl II2000I) study. The 50% completeness is 2 Mq for this sample, as shown by the arrow. 
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Fig. 9. — The derived IMF's between 3—7 pc assuming a 2 Myr isochrone (top), a 4 Myr isochronc (middle) and an average of the IMF's 
derived for the 2 to 4 Myr isochroncs in 0.5 Myr increments (bottom). The error bars on the top two panels include Poisson noise and the 
uncertainties in the completeness corrections. The number of stars in each mass bin in the average IMF is determined as the average of 
the 5 IMFs for ages 2-4 Myr. The error bars for each bin in the average IMF are then calculated as the standard deviation of the number 
of stars in the bin. The arrow shows the 50% completeness limit in each panel and the number next to arrow indicate the limiting mass 
in solar masses. For the average IMF, the 50% completeness limit is determined by the 4 Myr isochrone since here the completeness limit 
corresponds to the highest mass. 
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Fig. 10. — An illustration of the effect on the derived IMF of individual random differential extinction within 30 Dor. The solid line 
is an artificial Salpeter mass function and the dotted line is the retrieved mass function by adopting a magnitude-limited sample. The 
magnitude limit was assumed to be the same as used bv lSirianni et al.l I I2000I ) {V=2A.7 mag). Crosses show the mass function derived by 
[STrianni ct al. (2000). The stair-case shape of the input IMF at the high-mass end is due to finite, yet constant size of the mass bins. The 
50% completeness limit at 1.35 Mq is indicated by the arrow. 
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Fig. 11. — A comparison oftwo cumulative luminosity distributions outside 3 pc and down to the 50% completeness level of F160W= 21 
mag for the 3—5 pc annulus. A K— S test of the two distributions gives a maximum distance between the two distributions of 0.039 and the 
probability the two distributions to be drawn from the same parent distribution is 10%. The implication is that there is no evidence for 
luminosity segregation in the two annuli. 
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Fig. 12. — The surface brightness profile of the R136 cluster within a 7 pc radius of the cluster center. The solid line is a pure power— law 
fit from 0.02 pc to 2 pc with a derived slope of —1.48. The dashed line is an Elson, Fall & Freeman type profile with a core radius of 
0.025pc and a power— law slope of —1.54. The fit was done from O.OOQpc— 2pc. Beyond 2 pc, the presence of the individual bright stars 
labeled in Fig. |2] introduces the jitter seen in the surface brightness profile. 



The low-mass IMF in the 30 Doradus cluster 



23 



Annulus(pc) 


slope(Q) 


Max{F160W) 


0.6- 1.0 


0.38 ±0.07 


18.0 


1.0-2.0 


0.32 ± 0.03 


19.5 


2.0-3.5 


0.31 ± 0.02 


20.5 


3.5 - 5.0 


0.32 ± 0.01 


21.0 


5.0-7.0 


0.31 ±0.01 


21.5 



TABLE 1 

The derived slopes ( dFieow °^ 10"^^®°^) of the F160W band luminosity functions for R136 together with their 

UNCERTAINTIES. THE INNER AND OUTER RADII ARE GIVEN FOR EACH ANNULUS. THE MAXIMUM F160W BAND MAGNITUDE (50% 

completeness limit) USED TO DERIVE THE FIT IS SHOWN AS WELL. 



Xcen 


Ycen 


F160W 


FieOWerror 


104.546 


280.899 


19.09 


0.04 


156.263 


148.125 


20.28 


0.05 


121.580 


148.709 


20.15 


0.03 


315.194 


150.819 


20.29 


0.03 


481.127 


152.319 


20.27 


0.04 


69.890 


154.680 


17.73 


0.05 


280.115 


155.118 


19.76 


0.01 


91.993 


157.119 


18.87 


0.04 


443.582 


157.392 


17.34 


0.02 



TABLE 2 

The F160W band detected sources in R 136. The full table is available online. 



)(Myr) 


Annulus(pc) 


Mass ran; 


?e(Mp,) 


Slope(r) 


Mass range(M0) 


Slope(r) 


3 


0.6-1.0 


8.9- 


20 


-1.7 ±0.3 






3 


1.0-2.0 


8.9- 


20 


-1.5 ±0.1 






3 


2.0-3.0 


8.9- 


20 


-1.4 ±0.2 






3 


3.0-5.0 


1.4- 


20 


-1.2 ±0.1 






3 


5.0-7.0 


0.8- 


20 


-0.9 ±0.1 


1.4-1.7 


-0.9 ±0.2 


2 


3.0 - 5.0 


1.1 - 


20 


-1.0±0.1 


1.1-1.7 


-1.3 ±0.3 


2 


5.0 - 7.0 


0.7 - 


20 


-1.0 ±0.1 


0.7-1.7 


-0.8 ±0.2 


4 


3.0 - 5.0 


1.4 - 


20 


-1.2±0.1 






4 


5.0 - 7.0 


1.1 - 


20 


-1.2±0.1 


1.1-1.6 


-1.3 ±0.4 


- 4 


3.0 - 5.0 


1.4 - 


20 


-1.3±0.1 






- 4 


5.0 - 7.0 


1.1 - 


20 


-1.2±0.1 


1.1-1.6 


-1.1 ±0.4 



TABLE 3 

The derived slopes (r, dN/dlogM <x M^) of the mass functions together with the derived uncertainties of the slopes. The 
slope of a Salpeter IMF in these units is —1.35. All the fits were performed over the mass range indicated in the table. 



Center of mass bin (Mq) 
Annulus(pc) 


0.7 


0.9 


1.1 


1.4 


2.0 


2.8 


3.5 


5.0 


7.1 


8.9 


11.2 


14.1 


17.8 


22.4 


28.2 


0.6-1.0 














2 


11 


13 


21 


42 


51 


50 


25 


18 


11 


8 


1.0-2.0 


9 


27 


44 


56 


92 


155 


85 


56 


111 


91 


77 


48 


35 


30 


10 


2.0-3.0 


40 


74 


126 


139 


163 


157 


100 


53 


105 


80 


55 


57 


30 


17 


9 


3.0-5.0 


283 


328 


416 


380 


329 


268 


137 


86 


168 


111 


73 


41 


40 


29 


9 


5.0-7.0 


310 


351 


355 


297 


278 


202 


102 


56 


108 


74 


47 


43 


24 


24 


6 



TABLE 4 

The number of stars per mass bin used to derive the IMF for the 3 Myr isochrone. The numbers have not been corrected 

for incompleteness. 



